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CONVERSION  FACTORS,  NON- SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 


BY 


To  Obtain 


cubic  feet 
degrees  (angle) 

Fahrenheit  degrees 
feet 

pounds  (force) 

pounds  (force)  per  square  inch 
pounds  (mass)  per  cubic  foot 
seconds  (angle) 
seconds  (sidereal) 
square  feet 


0.02831685  cubic  metres 

0.01745329  radians 

5/9  (F  -  32)  Celsius  degrees  of  Kelvins1 

0.3048  metres 

4.448222  newtons 

0.006894757  megapascals 

16.01846  kilogram  per  cubic  metre 

0.000004848  radians 

0.9972696  seconds 

0.09290304  square  metres 


1  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  <F)  readings, 
use  the  following  formula:  C  -  <5/9)(F  -  32).  To  obtain  Kelvin  <K)  readings, 
use:  K  -  (5/9) (F  -  32)  •*>  273.15. 


EFFECTS  OF  EXPLOSIONS  IN  UNDERGROUND  MAGAZINES 


CHAPTER  I 
INTRODUCTION 

1.1  OVERVIEW 

The  advantages  of  underground  facilities  for  storage  of  explosives 
and  ammunition  include  constant  humidity  and  temperature  conditions,  and 
greater  protection  against  enemy  attack,  fire,  intrusion,  or  sabotage. 
However,  the  establishment  of  adequate  safety  standards  (assuming 
accidental  explosions  in  the  magazines)  must  precede  their  use.  These 
safety  standards  must  identify  the  required  depth  of  rock  or  soil  cover, 
minimum  spacing  between  storage  chambers,  and  the  minimum  safe  distance 
between  the  portal  and  inhabited  buildings  and/or  public  roads.  These 
parameters  are  dependent  upon  the  geologic  media,  underground  geometrical 
configurations,  type  and  quantity  of  explosive,  and  the  topography 
outside  the  portal.  Since  the  number  cf  variables  is  extensive,  a 
comprehensive  analytical  and  experimental  analysis  of  the  problem  is 
prohibitive.  As  a  result,  several  aspects  of  the  problem  had  not  been 
investigated  adequately  at  the  time  of  this  study. 

The  major  portion  of  the  work  reported  here  was  conducted  in  1979-80, 
and  the  "previous"  studies  referred  to  in  the  text  of  this  report  were 
published  prior  to  1979 r 

1.2  BACKGROUND 

In  most  instances,  the  greatest  hazard  from  accidental  explosions  of 
deeply  buried  magazines  is  airblast.  Consequently,  previous 
investigators  have  addressed  primarily  the  airblast  standoff  distance 
required  for  inhabited  buildings.  The  general  approach  has  been  to 
obtain  empirical  mathematical  relationships  based  upon  small-scale  model 
test  data  to  predict  the  peak  pressure  generated  at  the  portal  of  a 
particular  geometry  as  a  function  of  the  explosive  loading  density  of  the 
magazine.  Empirical  relationships  are  then  developed  to  predict  the  peak 
free -field  pressures  at  scaled  distances  based  upon  the  peak  pressure  at 
the  portal.  For  the  geometrical  configuration  and  range  of  loading 
densities  investigated,  this  predictive  technique  is  adequate;  however. 


extrapolation  of  the  empirical  relationships  to  different  geometries 
does  not  appear  to  be  mathematically  rigorous  and  may  not  produce  valid 
estimates.  These  investigations  will  be  discussed  in  more  detail  in 
Chapter  2. 

The  goal  is  to  suppress  the  explosion  and  to  reduce  the  standoff 
distance  required  for  buildings  as  a  result  of  the  ground  shock  or 
airblast  hazard.  However,  prior  to  this  investigation,  significant 
deficiencies  existed  in  prediction  of  ground  shock  near  decoupled2 
cylindrical  cavities  of  finite  length.  The  ground  shock  produced  by 
decoupled  cylindrical  charges  had  been  investigated  only  at  large 
distances  from  the  cavity.  The  effects  of  decoupling  near  a  cavity  had 
been  investigated  only  for  spherical  geometries.  As  a  result,  the 
separation  distances  between  storage  chambers  required  to  prevent 
explosive  communication  could  not  be  determined  accurately.  Similarly, 
it  was  not  possible  to  predict  the  chamber  separation  required  to  prevent 
an  explosion  in  one  chamber  from  damaging  the  concents  of  an  adjacent 
chamber.  Ideally,  only  minor  surface  damage  to  the  wall  of  the  nearest 
adjacent  chamber  should  occur.  And  finally,  additional  Information  was 
needed  in  order  to  determine  the  safe  standoff  distance  between 
underground  magazines  and  residential  buildings  required  to  preclude 
structural  damage  as  a  result  of  ground  motion. 

1.3  OBJECTIVE 

The  general  objective  of  this  investigation  was  to  improve  the 
state- of >thc» art  of  prediction  techniques  used  for  establishing 
underground  explosive  storage  safe‘s  standards.  Since  the  ground  shock 
produced  by  decoupled  explosions  in  cylindrical  cavities  was  not  well 
defined,  emphasis  was  on  determining  standoff  distances  required  to 
prevent  explosive  communication,  separation  required  to  minimize  damage 
to  contents  in  adjacent  chambers,  and  safe  standoff  distances  for  ground 
shock  effects  on  public  buildings  and  highways.  A  secondary  objective 
was  to  determine  safe  standoff  distances  for  airblast  propagating  from 
the  postal. 


2  A  "decoupled"  charge  Is  defined  as  a  cavity  or  chamber  loaded  with 
explosives  which  occupy  less  than  the  entire  volume  of  the  chamber. 
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1.4  APPROACH 


The  following  approach  was  used  to  develop  equations  to  predict  safe 
standoff  distances.  First,  literature  germane  to  each  major  aspect  of 
the  problem  was  reviewed.  The  gas  dynamics  part  of  the  problem  included 
the  peak  pressure  generated  in  vented  and  unvented  cavities,  blast 
propagation  through  tunnels  and  free-field  blast  propagation.  The  ground 
motion  survey  included  analytical  and  experimental  investigations  of  the 
strain  and  velocity  fields  produced  by  decoupled  explosions  in  cavities 
in  rock,  and  the  analysis  of  the  spall  process  at  free  rock  surfaces. 

The  literature  survey  is  presented  in  Chapter  2. 

Second,  experimental  data  voids  were  identified  which  precluded  an 
analysis  of  the  cylindrical  decoupling  process.  A  small-scale  model  was 
designed  and  used  to  obtain  the  additional  data  required  to  predict  the 
freo-field  strain  and  velocity  field.  Additionally,  the  model  was 
designed  to  measure  Che  free-field  airblast  pressures  generated  by  the 
blast  exiting  the  portal.  The  design  of  the  model  is  discussed  in 
Chapter  3, 

Third,  the  data  were  analysed  and  the  ground  shock  propagation 
process  was  deduced  for  decoupled  detonations  in  cylindrical  cavities  of 
finite  length  (Chapter  4).  Equations  were  developed  which  express  the 
relationship  between  free-field  strain,  the  geometry  of  the  cavity,  the 
properties  of  the  media  and  the  loading  density  of  the  chamber. 

Fourth,  equations  v ere  developed  to  predict  the  respective  ground 
shock  standoff  distances  required.  Those  equations  were  compared  to 
those  currently  employed  in  the  1980  Department  of  Defense  (D00) 
Ammunition  and  Explosive  Safety  Standards  (Reference  1).  The  details  are 
presented  in  Chapter  5. 


CHAPTER  2 
LITERATURE  REVIEW 


Two  major  explosion  phenomena- -gas  dynamics  and  ground  shock--are 
relevant  to  the  determination  of  standoff  distances  required  for 
underground  magazines.  The  gas  dynamics  area  includes  the  detonation 
pressure  generated  in  the  storage  chamber,  attenuation  of  the  shock  front 
during  propagation  through  the  magazine  complex,  prediction  of  the  blast 
pressure  history  at  the  portal,  and  the  external,  free-field  decay  of  the 
blast  wave.  The  ground  shock  area  includes  determination  of  the:  stress 
induced  in  the  storage  chamber  walls,  analysis  of  the  free-field,  and 
determination  of  the  response  of  free  surfaces  as  a  function  of  the 
incident  strain.  Relevant  literature  in  each  of  these  areas  will  be 
discussed. 

2.1  GAS  DYNAMICS 

2.1.1  Cavity  Pressure.  Explosions  produce  instantaneous  release  of 
tremendous  amounts  of  energy.  Swung  shock  waves  quickly  develop  and 
interact  with  the  cavity  boundaries  (Brode  and  Parkin,  Reference  2). 

After  numerous  reflections  and  shock  interactions,  the  waves  are  rapidly 
dissipated  and  a  quasistatic  pressure  is  generated,  referred  to  here  as 
the  chamber  pressure  (Proctor  1972,  Reference  3).  The  peak  wall 
pressures  developed  by  reflection  of  the  initial  shock  can  be 
significantly  (8  to  10  times)  larger  than  the  chamber  pressure 
(Skjeltorp,  Hegdahl,  and  Janssen,  Reference  4).  Although  the  reflected 
shocks  have  relatively  sraal*  impulse,  localized  cracking  of  the  Cumber 
and  nearby  tunnel  walls  ca  i  be  produced.  In  a  closed  u^ivity,  the  chamber 
pressure  decays  slowly  as  a  result  of  heat  loss  to  the  walls.  If  the 
cavity  is  vented,  the  chamber  pressure  decay  rate  is  controlled  by  the 
magnitude  of  the  mass  flux  of  the  detonation  products  (Proctor, 

Reference  5). 

Filler  (Reference  6)  developed  a  formulation  to  predict  the  chamber 
pressure,  using  an  elementary  energy  balance  which  assumed  that  the  ratify 
of  specific  heat  was  constant  and  that  an  oxygen- deficient  condition  did 
not  exist.  Consequently,  this  analysis  was  limited  to  low  loading 
densities.  For  TNT,  the  technique  is  valid  fur  chamber  loa^'"^  densities 


less  then  0.02  lb/ft3.  An  empirical  fit  to  experimental  data  for  TNT  was 
suggested  by  Weibul  (Reference  7).  However,  the  equation  was  limited  to 
the  range  of  the  reported  data,  and  there  was  not  an  effective  means  of 
relating  TNT  data  to  other  explosives. 

Theoretical  techniques  have  also  been  developed  to  predict  the 
chamber  pressures  generated  by  explosions  in  cavities.  Proctor  (1974, 
Reference  5)  wrote  a  flexible  computer  code  for  explosives  containing 
carbon,  hydrogen,  oxygen  and/or  aluminum.  Based  upon  the  chemical 
composition  of  the  explosive  and  the  quantity  of  air  in  the  chamber,  a 
chemical  balance  is  used  to  determine  the  quantities  of  Al2,  03,  H20,  CO, 
C02,  C,  H2,  02,  and  N2  that  are  generated  by  the  chemical  reaction.  The 
amount  of  thermal  energy  (Q)  generated  by  the  reaction  was  then 
calculated.  To  facilitate  computations,  Proctor  assumed  that  the  mixture 
of  detonation  products  was  initially  at  ambient  temperature  and  pressure. 
The  thermal  energy  released  by  the  reaction  is  added  to  the  ambient 
condition  by  a  constant  volume  process, 

Q  -  (T)  dT 

where  Cv  is  the  specific  heat  of  the  gas  at  constant  volume  and  T  is  the 
equilibrium  temperature . 

The  integral  equation  is  solved  by  iteration  for  the  equilibrium 
temperature  of  the  final  mixture  which  satisfies  the  equation.  Using  the 
equilibrium  temperature,  mass  of  products,  and  the  volume  of  the  chamber, 
an  equation  of  state  was  employed  to  calculate  the  chamber  pressure. 
Figures  2.1  and  2.2  provide  comparisons  of  Proctor's  (1974,  Reference  5) 
theoretical  predictions  with  experimental  data  for  TNT  and  PETN 
explosives.  Stromsoe  (Reference  8)  calculated  the  temperature  and 
chamber  pressure  developed  in  an  unvented  cavity  by  the  detonation  of  TNT 
charges.  For  the  range  of  loading  densities  investigated  (0.012  to  16,88 
lb/ft3),  his  results  were  identical  to  Proctor's, 

2.1.2  Blast  Propagation  in  Tunnels.  Following  an  explosion  in  an 
underground  magazine ,  venting  of  the  detonation  products  through  the 
chamber  entrance  results  in  the  propagation  of  shock  waves  through  the 
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access  tunnel.  The  physical  phenomena  involved  are  complex,  and  a 
generalized  prediction  of  the  blast  pressures  produced  inside  the  tunnel 
or  at  the  portal  is  very  difficult. 

First,  consider  the  movement  of  the  blast  wave  into  the  access  tunnel 
from  the  storage  chamber.  Upon  entering  the  tunnel,  the  peak 
overpressure  is  significantly  reduced  by  diffraction  (Porzel, 

Reference  9).  As  the  shock  front  turns  and  the  flow  is  accelerated,  a 
pressure  drop  occurs.  Diffraction  is  quickly  followed  by  other 
irreversible  losses.  A  few  diameters  inside  the  tunnel,  the  peak  shock 
pressure  increases  over  the  diffraction  pressure.  This  basically  results 
from  the  interrelated  process  of  preferential  flow,  multiple  reflections 
and  shocking  up.  Deeper  in  the  tunnel  (within  approximately  20 
diameters),  the  turbulent  boundary  layers  of  the  shock  front  converge, 
and  a  turbulent  choke  develops.  The  turbulent  choke  is  believed  to  form 
at  a  constant  distance  behind  the  shock  front.  The  flow  then  becomes  a 
quasi-steady  state  process. 

The  roughness  of  the  tunnel  wall  surface  also  attenuates  the  peak 
blast  pressures  (Skjeltorp  1968,  Reference  10).  For  a  typical 
underground  storage  facility,  surface  roughness  could  reduce  the  pressure 
outside  the  portal  by  a  factor  of  two  (or  more),  compared  to  an 
essentially  smooth-walled  facility  (Skjeltorp  1975,  Reference  11).  The 
degree  of  attenuation  depends  upon  the  shock  strength,  as  well  as  the 
dogree  of  surface  roughness  (Kriebel,  Reference  12). 

Porzel *s  theoretical  analysis  and  several  experimental  investigations 
have  demonstrated  that  macroscopic  surface  roughness  in  tunnels 
significantly  attenuates  weak  shocks  by  surface  friction.  In  contrast, 
strong  shocks  are  relatively  insensitive  to  macroscopic  surface 
roughness.  Even  if  there  are  no  viscous  losses  to  surface  friction, 
attenuation  in  straight  tunnels  still  occurs  as  a  result  of  rarefaction. 
Advancement  of  the  rarefaction  wave  toward  the  shock  front  results  in  a 
reduction  of  the  peak  pressure  and  an  increase  in  the  duration  of  the 
positive  phase.  Data  from  modal  tests  have  quantified  this  effect  as  a 
function  of  tunnel  length  (Schmidt,  Reference  13). 
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Many  underground  storage  complexes  are  much  more  complicated  than  a 
single  storage  chamber  connected  to  a  straight  access  tunnel.  Each  bend, 
corner,  or  intersection  in  a  tunnel  constitutes  a  geometrical 
discontinuity  which  impedes  shock  propagation  by  producing  multiple  shock 
reflections.  In  general,  the  more  complex  the  geometry,  the  greater  the 
attenuation  of  the  blast  pressure  as  it  propagates  through  a  tunnel. 

Theoretically,  blast  pressures  generated  in  an  underground  tunnel 
complex  can  be  described  mathematically  and  predicted  numerically 
(Skjeltorp,  Hegdahl,  and  Jenssen,  1975,  Reference  4).  In  most  studies  to 
date,  however,  the  approach  has  been  to  determine  experimentally  the 
blast  pressures  generated  in  small-scale  models  of  the  particular 
geometry  of  interest.  Semi-empirical  equations  are  fit  to  the  measured 
peak,  side-on  pressure  data  obtained  for  a  particular  geometry,  such  as 
that  shown  in  Figure  2,3.  According  to  Skjeltorp  (1968,  Reference  10), 
for  a  storage  chamber  connected  to  a  single  access  tunnel,  the 
relationship  between  peak  pressure,  charge  weight  in  the  magazine,  and 
dimensions  of  the  complex  is 

?  ■ 953  (C(^r,o-iis(^)s  0,45 

where  P  is  peak  side-on  pressure  (psi) 

Q  is  total  equivalent  TNT  charge  weight  (lbs) 

Aj  is  the  cross-sectional  area  of  the  access  tunnel  (ft2) 

Ag  is  tho  area  of  the  storage  chamber  (ft2) 
and  Vt  is  the  total  volume  (ft9)  of  the  chamber  and  tunnel  segment  to 
the  point  of  interest 

Peak  pressure  attenuation  rates  have  been  investigated  for  several 
basic  tunnel  dosigns  (Skjeltorp,  et  al,  1975,  Reference  4).  The  lumped 
approach  shown  above  illustrates  the  approach  used  in  most  experimental 
studies.  The  loading  densities  of  the  storage  chambers  investigated  have 
boon  less  than  6.3  lb/ft9,  with  most  less  than  3.25  lb/ft9. 

This  method  produces  very  good  results  for  the  geometry  investigated 
but  has  several  limitations,  according  to  Skjeltorp  (Reference  4).  The 
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data  obtained  is  applicable  only  to  geometrically  similar  models  (or 
full-scale  designs),  and  only  for  the  range  of  loading  densities  that 
were  investigated. 

2.1.3  Free-Field  Blast  Pressure.  An  accidental  explosion  in  an 
underground  magazine  normally  produces  a  hazardous  blast  wave  that 
propagates  outward  from  the  tunnel  exit.  This  external  blast  wave 
attenuates  at  different  rates  according  to  direction  from  the  exit;  i.e., 
as  a  function  of  angular  departure  from  the  extended  tunnel  center  line. 
Generally,  hazardous  airblast  effects  will  extend  much  farther  from  the 
exit  portal  than  will  the  detrimental  effects  of  ground  shock  or  debris. 
The  minimum  safe  distance  (or  airblast  standoff  distance)  is  usually 
defined  as  the  distance  from  the  portal  to  a  point  at  which  the 
overpressure  is  equal  to  0.725  psi  (Skjeltorp,  Reference  4). 

The  airblast  standoff  distance,  for  a  given  pressure  criterion,  is 
usually  predicted  from  empirical  relations  which  describe  the  attenuation 
of  the  peak  pressure  exiting  the  tunnel  portal;  i.e.,  the  exit  pressure. 
Fredrickson  and  Jenssen  (Reference  14)  obtained  free -field  airblast  data 
from  small-scale  model  tests  where  the  access  tunnel  exited  onto  a 
30 -degree  (downward)  slope.  Based  upon  these  results,  they  proposed  the 
empirical  relationship 

dc  -  k0Po0-67  (2.2) 

where  d5  is  the  airblast  standoff  distance 

P0  is  the  overpressure  at  the  tunnel  exit 
kg  is  a  constant  determined  by  the  geometrical  design  of  the 
complex  and  direction  of  interest  outside  the  tunnel. 

Fredrickson  and  Jenssen  noted  that  the  technique  may  not  yield 
accurate  quantity  distances  for  other  geometrical  configurations. 
Skjeltorp  (1975,  Reference  4)  subsequently  reviewed  Fredrickson  and 
Janssen's  data  and  proposed  the  empirical  equation 

-  CePo0-57  «-3) 
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where  r  is  the  radial  distance  in  feet  from  the  portal,  along  a  radial 
line  originating  at  the  portal  to  a  point  at  which  the  peak 
blast  pressure  has  decayed  to  0.725  psi. 

D5  is  the  diameter  of  the  access  tunnel 

Cq  is  a  directional  coefficient  functionally  dependent  upon  the 
direction  of  measurement  with  respect  to  the  extended  centerline 
and  6  is  the  angle  between  the  direction  of  the  extended  centerline. 

The  values  of  C9  as  a  function  of  6  are: 

6(deg) _ 0 _ 30 _ 60 _ 90  120 

C9  3.00  2.67  2.00  1.33  0.75 

Skjeltorp  suggests  that  Equation  2.3  is  valid  for  any  underground 
geometry  if  the  exit  geometry  is  equivalent  to  that  used  by  Fredrickson. 

In  a  study  involving  a  dual  storage  chamber  design,  Gurke  and 
Scheklinski  (Reference  15)  obtained  free -field  blast  data  along  the 
extended  centerline  of  the  main  access  tunnel.  The  model  was  fabricated 
from  steel  tubes  which  exited  onto  a  flat  surface  (0-degree  exit 
surface).  At  low  loading  densities,  the  peak  blast  pressures  were 
substantially  lower  than  those  predicted  by  Equation  2.3.  At  larger 
loading  densities  (1.0  to  2.5  lb/ft3),  the  blast  data  obtained  along  the 
extended  centerline  access  agreed  with  Equation  2.3.  Thus ,  the  rate  of 
decay  along  the  extended  centerline  was  not  affected  by  exit  slopes 
between  0  and  30  degrees.  These  results  provided  no  additional 
information  about  the  directional  sensitivity  of  Equation  2.3  to  the  exit 
geometry.  The  blast  pressures  generated  along  the  extended  centerline 
should  be  less  sensitive  to  the  exit  geometry  than  any  other  direction. 

Several  investigations  have  been  conducted  to  determine  the  external, 
free-field  blast  pressures  generated  by  gun  firings.  Westine 
(Reference  16)  reviewed  data  obtained  over  a  wide  range  of  operating 
conditions,  from  fourteen  different  weapons.  The  barrel  lengths  and  bore 
diameters  varied  more  than  an  order  of  magnitude.  The  total  energy  of 
the  propelling  charges  varied  by  a  factor  of  4490.  Peak  free-field  blast 
data  were  obtained  along  radial  lines  in  a  horizontal  plane  that  was  an 
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extension  of  the  bore.  Each  of  the  measurements  made  along  30  to 
90-degree  lines  showed  the  peak  blast  pressure  decaying  approximately 
as  r‘3/2. 

A  comparison  of  the  results  obtained  by  the  various  investigators 
(Skjeltorp,  Fredricksen,  Jenssen,  Gurke,  Scheklinski,  and  Westine) 
suggests  that  the  directional  coefficients  in  Equation  2.3  in  the  30,  60 
and  90-degree  directions  are  not  sensitive  to  differences  in  exit  slopes 
between  0  and  30  degrees. 

2.2  GROUND  SHOCK 

It  has  been  shown  both  theoretically  (Herbst,  Werth,  and  Springer, 
Reference  17)  and  experimentally  (Perret,  Reference  18)  that  decoupling 
the  explosive  energy  by  detonating  charges  in  cavities  larger  than  the 
charge  volume  is  an  effective  method  of  reducing  ground  shock  magnitudes. 
The  decoupling  in  spherical  cavities,  for  example,  has  been  explored  as  a 
possible  method  for  concealing  underground  nuclear  explosions.  The 
technique  is  also  used  in  the  blasting  industry  in  presplitting,  smooth 
wall,  and  cushion  blasting  operations. 

Since  munitions  and  explosives  typically  occupy  only  a  few  percent  of 
the  total  volume  of  underground  magazines,  the  net  explosive  weight  of  an 
accidental  explosion  would  be  comparable  to  a  concentrated  charge  that  is 
highly  decoupled.  Thus,  the  amount  of  decoupling  is  inversely 
proportional  to  the  magazine  loading  density. 

Decoupling  should  result  in  much  less  energy  coupled  into  the  earth 
by  accidental  explosions  in  magazines  than  would  be  predicted  based  on 
data  from  fully-coupled  tests.  Theoretical  investigations  of  decoupling 
in  cylindrical  cavities  have  not  been  published.  However,  empirical  data 
(Atchison,  Duvall,  and  Pugliese,  Reference  19)  were  obtained  for  radial 
distances  greater  than  10  cavity  lengths,  for  matching  decoupled  and 
fully-coupled  cylindrical  charge  detonations.  Drake  (Reference  20) 
demonstrated  that,  at  these  ranges,  the  reduction  in  ground  shock 
achieved  by  decoupling  cylindrical  charges  is  roughly  equivalent  to  that 
of  spherical  charge  decoupling.  The  peak  particle  velocity  obtained  for 
decoupled  charges  (v)  divided  by  the  peak  particle  velocity  obtained  for 
fully  coupled  charges  (vQ) ,  expressed  as  a  function  of  the  relative 
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loading  density,  is  shown  in  Figure  2.4.  The  data  for  granite  and 
limestone  were  obtained  using  cylindrical  cavities  containing  chemical 
explosives.  Data  for  spherical  cavities  in  salt  involved  both  chemical 
and  nuclear  explosives. 

The  velocity  and  strain  fields  produced  near  fully-coupled 
cylindrical  charges  have  not  been  reported.  A  small  quantity  of  close-in 
data  from  decoupled  cylindrical  charges  is  available  from  the  CIST 
program  (Amend,  References  21  and  22).  In  this  program,  a  5-lb/ft  linear 
charge  was  detonated  in  2- ft  diameter  vertical  holes  drilled  through 
layered  media.  Peak  particle  velocities  were  obtained  in  each  layer  at 
radial  distances  of  1.5  to  6  cavity  diameters.  Five  peak  particle 
velocity  measurements  were  obtained  in  sandstone,  and  six  in  granite. 

The  quantitative  behavior  of  the  peak  strain  (or  particle  velocity) 
around  the  detonation  of  linear  charges  in  rock  can  be  deduced.  Near  a 
linear  charge  of  length  L,  the  stress  (or  strain)  field  along  the  length 
of  the  charge  is  determined  by  cylindrical  wave  propagation.  As  the 
cylindrical  wave  advances,  the  amplitude  of  the  wave  front  decays  as 
r"1/2,  where  r  is  the  radial  distance  from  the  center  of  the  charge 
(Person,  Reference  23).  At  large  distances  from  a  linear  charge  (i.e., 
where  r»L) ,  the  stress  field  approximates  that  produced  by  a  spherical 
charge.  Ideally,  the  amplitude  of  spherical  waves  decays  as  r*1. 

However,  as  a  result  of  irreversible  losses  (heat,  mechanical  crushing, 
etc.),  the  decay  in  rock  is  more  nearly  proportional  to  r'a  (Cooper, 
Reference  24).  As  indicated  in  Figure  2.5,  between  those  limits,  a 
transition  from  cylindrical  to  spherical  decay  must  occur. 

2.2.1  Tunnel  Damage.  Underground  installations  con  be  significantly 
damaged  by  the  stress  waves  generated  from  explosions  within  the 
facility.  The  magnitude  of  damage  is  dependent  on  the  energy  released  by 
the  charge,  the  degree  of  charge  coupling  (loading  density),  properties 
of  the  medium,  and  the  distance  between  the  charge  and  adjacent  tunnels. 
Several  investigations  have  been  conducted  to  determine  the  strain  (or 
velocity)  threshold  at  which  various  types  of  tunnel  damage  occurs. 

Damage  to  tunnels  is  generally  divided  into  four  categories,  based  upon 
the  damage  profile  (see  Figure  2.6)  along  the  tunnel  wall  nearest  the 
charge.  In  Army  nuclear  weapons  effects  manuals  (e.g.,  Reference  27), 
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the  four  zones  (1,  2,  3  and  4)  are  called  the  zone  of  complete  damage, 
the  zone  of  rock  breakage,  the  zone  of  continuous  slabbing,  and  the  zone 
of  discontinuous  damage,  respectively. 

The  Underground  Test  (UGT)  Program  (Engineering  Research  Associates, 
References  28  and  29)  provided  damage  data  in  three  types  of  rock 
(limestone,  granite  and  sandstone).  The  tunnels,  in  jointed  granite  and 
sandstone,  had  horseshoe -shaped  cross-sections.  The  diameters  of  the 
tunnels  tested  were  6,  15,  and  30  ft.  Fully-coupled  TNT  charges,  varying 
in  weight  from  320  to  32,000  pounds,  were  detonated  at  various  distances 
from  the  test  tunnel. 

In  addition,  a  series  of  153  model  tests  were  conducted  at  the 
Waterways  Experiment  Station  in  rock-matching  grout  having  nominal 
unconfined  compressive  strengths  of  10,000,  4,000  and  500  psi  (Joachim, 
Reference  30) .  Various  standoff  distances  and  tunnel  diameters  were  also 
tested  in  these  models,  Hendron,  Clark,  and  Strange  (Reference  31)  used 
a  model  which  had  eight  tunnels  of  various  diameters  oriented  about  a 
2 -lb  TNT  spherical  charge.  The  distances  from  the  center  of  the  charge 
to  the  centerline  of  the  tunnels  ranged  from  1  to  5  ft.  Other  data  were 
obtained  in  a  series  of  23  tests  conducted  In  Kayenta  sandstone  near 
Grand  Junction,  CO,  Core  samples  indicated  a  relatively  uniform  rock 
with  thin,  horizontal,  irregularly  spaced,  clay  seams  (Swift, 

Reference  32).  Stemmed  charges  were  detonated  at  various  distances 
from  2,  3,  and  4* ft  diameter  drill  holes. 

Based  upon  the  results  of  these  Investigations,  the  following 
conclusions  are  implied:  Zone  1  failure  occurs  when  the  compressive 
stress  exceeds  the  ultimate  compressive  strength  of  the  partially- 
confined  medium  adjacent  to  the  opening;  this  usually  corresponds  to  a 
strain  of  0.5  to  1.0  percent.  Zone  2  failure  occurs  when  the  local 
compressive  strain  around  the  opening  equals  the  strain  at  ultimate 
conditions  in  partially-confined  compression.  For  most  rocks,  this 
corresponds  to  a  strain  between  0.2  and  0.4  percent.  Zone  3  damage 
occurs  when  the  local  compressive  strain  around  the  opening  exceeds  the 
strain  corresponding  to  the  ultimate  strength  in  unconfined  compression. 
Typically,  this  strain  is  between  0.06  and  0.2  percent.  Zone  4  damage 
develops  when  the  local  tensile  strain  around  the  opening  equals  the 
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strain  corresponding  to  the  effective  tensile  strength  of  the  rock. 

For  most  types  of  rock  this  type  damage  would  occur  at  free-field  strains 
of  0.03  to  0.06  percent. 

2.2.2  Snail.  The  conditions  under  which  spalling  occurs  are  well 
understood  (Rinehart,  Reference  33).  When  a  compressional  wave 
intersects  a  free  surface,  it  will  be  reflected  as  a  tension  wave,  and 
the  total  stress  at  a  point  in  the  rock  wall  will  be  the  sum  of  the 
incident  compressioi  il  stress  and  the  reflected  tension  stress.  Whenever 
the  tension  exceeds  the  compression  by  an  amount  approximately  equal  to 
the  ultimate  tensile  strength  of  the  rock,  it  will  fracture  and  a  spall 
will  develop.  The  velocity  of  the  spall  will  depend  upon  the  magnitude 
of  the  momentum  trapped  within  the  spall. 

If  the  net  stress  level  is  greater  than  twice  the  ultimate  tensile 
strength  of  the  rock,  multiple  spalls  will  develop.  Assume  that  the 
first  spall  was  created  as  indicated  above.  The  remaining  part  of  the 
compressional  wave  will  strike  the  new  free  surface,  and  the  process  will 
be  repeated  until  the  amplitude  of  the  tensile  stress  is  less  than  the 
ultimate  tensile  strength  of  the  rock. 

Fractures  in  the  medium  may  significantly  affect  spalling 
(Reference  33).  If  the  fractured  surfaces  are  in  intimate  contact;  a 
compression  wave  can  be  transmitted  through  the  interface  with  little  or 
no  degradation;  however,  the  interface  cannot  support  tension.  When  the 
reflected  tensile  wave  from  a  free  surface  returns  to  the  fracture,  a 
spall  will  occur  at  the  fracture.  The  velocity  of  the  spall  will  depend 
upon  the  momentum  trapped  within  the  spalled  slab  of  rock,  which  is 
dependent  on  the  slab  thickness. 

Spalling  in  real  materials  is  often  complicated  by  attenuation, 
divergence,  inhomogeneities  in  the  material  and  finite  fracture  times. 
Spall  velocity  predictions  are  often  based  upon  idealised  assumptions. 

In  a  safety  analysis,  this  is  not  practical  because  the  conditions  at  the 
free  surface  are  not  always  known,  and  upper  bound  predictions  may  not 
result.  It  appears  more  realistic  to  specify  the  required  separation  of 
underground  storage  chambers  in  terms  of  the  maximum  possible  spall 
velocity,  which  is  always  equal  to  or  less  than  twice  the  particle 
velocity  at  the  free  surface. 


20 


2.3  SUMMARY  (Literature  Review) 

The  interrelated  processes  that  govern  the  effects  of  explosions  in 
underground  magazines  have  not  been  sufficiently  investigated  to 
establish  comprehensive  safety  standards.  Theoretical  techniques  have 
been  developed  to  accurately  predict  the  pressure  generated  by  explosions 
in  underground  chambers.  Peak  blast  pressures  produced  in  the  connecting 
passageways  and  at  the  portal  can  be  predicted  empirically,  but  only  for 
the  limited  number  of  geometries  and  loading  densities  that  have  been 
tested. 

The  ground  shock  produced  by  decoupled  explosions  in  cylindrical 
cavities  also  has  not  been  adequately  investigated.  It  has  been 
demonstrated  that  decoupling  will  significantly  reduce  the  peak  strain 
and  velocity  field.  Near  cylindrical  cavities  of  finite  length,  where 
(r«L) ,  the  decay  in  the  peak  values  near  the  cavity  will  be  dictated  by 
cylindrical  wave  expansion,  Far  from  the  cavity,  where  (r»L),  spherical 
wave  expansion  will  predominate.  The  radial  distance  at  which 
cylindrical  decay  changes  to  spherical  has  not  been  determined 
experimentally.  Consequently,  the  close-in  stress  and  strain  fields 
cannot  be  predicted  with  any  degree  of  certainty. 

The  magnitude  of  damage  to  unlined  tunnels  as  a  result  of  the  stress 
field  generated  by  an  external  explosion  can  be  related  to  the  strain 
developed  at  the  tunnol  wall.  In  unlined  tunnels,  catastrophic  failure 
(tunnel  closure)  may  occur  if  the  strain  at  the  tunnel  wall  is  greater 
than  0.4  percent.  Catastrophic  wall  failure  of  a  magazine  could 
conceivably  result  in  explosive  initiation  of  the  contents.  Minor 
surface  damage  may  occur  at  strain  magnitudes  greater  than  0.03  percent. 
To  preclude  damage  to  the  contents  of  adjacent  underground  storage 
chambers,  the  chamber  separation  should  be  largo  enough  to  prevent  any 
wall  surface  damage  in  the  chambers  adjacent  to  an  explosion. 

From  the  literature  survey,  it  was  apparent  that  a  major  problem 
exists  in  the  current  ability  to  predict  the  magazine  spacings  required 
to  prevent  interchamber  damage  by  ground  shock.  No  close-in  ground  shock 
data  are  available  for  decoupled,  cylindrical  cavities.  A  small-scale 
model  test  program  was  therefore  designed  to  obtain  such  data  as  a  part 
of  this  study. 
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Figure  2,3.  An  example  of  the  scaling  procedure  employes  to  predict 
peak  pressures  in  underground  magaxinea* 


U 


<0 


X 

< 

111 


Figure  2.5.  Qualitative  behavior  of  the  strain  field  produced 
around  line  charges  of  finite  length. 
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Figure  2.6.  Typical  tunnel  damage  profile. 


CHAPTER  3 


EXPERIMENTAL  MODEL 

3.1  INTRODUCTION 

A  small-scale  model  (Figure  3.1)  was  used  to  measure  the  strain, 
velocity  and  acceleration  field  produced  in  rock  by  uncoupled  cylindrical 
explosions,  and  to  measure  the  free- field  blast  produced  outside  the 
portal.  From  this  data,  the  safe  standoff  distances  for  ground  shock  and 
air  blast  effects  could  be  inferred  for  accidental  detonations  in 
full-scale  magazines.  No  attempt  was  made  to  duplicate  a  particular 
full-scale  facility. 

3.2  DESCRIPTION  OF  MODEL 

The  model  consisted  of  a  single  storage  chamber  and  access  tunnel 
embedded  in  an  instrumented,  granite -matching  grout  block.  The 
dimensions  (diameter  and  length)  of  the  storage  chamber  correspond  to  a 
1: 75-scale  model  of  a  typical  underground  magazine.  A  comparison  of  the 
model  and  full-scale  dimensions  is  given  in  Table  3.1.  Surrounding  the 
storage  chamber  in  radial  arrays  were  six  accelerometers  and  eight  strain 
gages  (Figure  3.2).  The  distance  between  the  center  of  the  cavity  and 
each  ground  motion  gage  is  indicated  in  Table  3.2.  The  access  tunnel  was 
a  smooth-wall,  steel  tube  that  exited  the  block  tangential  to  an 
instrumented  concrete  slab.  The  scaled  diameter  of  the  access  tunnel  was 
chosen  to  match  the  diameter  of  a  typical  full-scale  design.  However, 
the  scaled  length  of  the  access  tunnel  was  somewhat  shorter  than  those 
used  in  typical  full-scale  designs.  Airblast  pressure  gages  were  placed 
(Figure  3.3)  along  four  lines  which  extended  in  different  angular 
directions  from  the  access  tunnel  portal.  Distances  between  blast  gages 
and  the  tunnel  portal  are  given  in  Table  3.3. 

3.3  DESIGN  CONCEPTS 

The  experiment  was  designed  to  measure  ground  shock  within  the  grout 
block  and  the  external  blast  pressures.  Because  of  the  absence  of 
existing  close-in  ground  shock  data  for  cylindrical,  decoupled 
explosions,  greater  emphasis  was  placed  on  the  acquisition  of  the  ground 
shock  data.  This  significantly  impacted  the  experimental  design. 
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The  relatively  simple  geometry  of  a  single  storage  chamber  and  access 
tunnel  was  used  for  two  reasons.  Multiple  storage  chambers  would  have 
increased  construction  cost  and  complicated  the  experiments.  The  single 
storage  chamber  design  also  eliminated  the  possibility  of  wave 
reflections  from  a  secondary  storage  chamber  that  could  have  compromised 
the  basic  ground  shock  measurements. 

The  length  of  the  storage  chamber  was  made  as  large  as  possible, 
relative  to  the  external  dimensions  of  the  model,  to  maximize  the  radial 
distance  dominated  by  the  cylindrical  stress  wave.  This  resulted  in  a 
relatively  short  access  tunnel,  compared  to  full-scale  designs. 

According  to  previous  studies,  however,  a  shorter  access  tunnel  would  not 
significantly  affect  the  usefulness  of  the  blast  pressure  data 
(Reference  4) . 

Reflections  from  the  outer  walls  of  the  model  limited  the  time 
duration  over  which  useful  ground  motion  data  could  be  obtained.  To 
minimize  this  problem  and  to  increase  the  radial  distance  available  for 
ground  shock  measurements,  the  storage  chamber  was  off-set  horizontally 
from  the  center  of  the  grout  block  (Figure  3.2).  The  time  required  for 
the  first  reflection  to  reach  each  sensor  was  calculated,  based  oh  the 
anticipated  grout  wave  speed.  These  results  indicated  that  the  transit 
times  were  sufficiently  large  that  the  measurement  of  peak  strain  or 
velocity  at  the  stations  selected  would  not  be  compromised.  The  accuracy 
of  the  predicted  reflections  are  discussed  in  Chapter  4. 

3.4  CONSTRUCTION 

After  the  exterior  concrete  forms  for  the  grout  block  were  set,  the 
forms  for  the  storage  chamber  and  access  tunnel  and  the  reinforcement 
material  were  suspended  by  wire.  Reinforcement  (Figure  3.4a)  was  used  to 
prevent  stress  concentrations  at  the  ends  of  the  storage  chamber  from 
cracking  the  block.  To  prevent  wave  reflections  from  compromising  the 
ground  motion  data,  very  little  reinforcement  was  used  in  the  region 
between  the  gages  and  the  storage  chamber  (Figure  3.4b). 

Since  thermodynamic  calculations  indicated  that  massive,  internal 
thermal  cracks  would  develop  if  the  test  block  was  formed  by  a  single 
pour,  three  equal  pours  were  made  on  three  consecutive  days.  Samples 
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were  taken  during  each  pour  for  later  laboratory  testing  to  confirm  the 
properties  of  the  grout.  The  sequence  of  the  pours  is  shown  in 
Figure  3.4.  Prior  to  the  second  and  third  pours,  the  previous  pour  was 
coated  with  epoxy  to  minimize  the  effects  of  cold  joints.  During  the 
second  pour,  thin-wall  copper  tubes  (Figure  3.4c)  were  inserted  for  later 
gage  placement.  To  insure  bonding  of  the  strain  gage  assembly  to  the 
test  block,  the  vertical  tube  over  the  center  of  the  storage  chamber  was 
removed  after  the  third  pour  reached  an  initial  set.  The  vertical  holes 
were  filled  with  same  grout  mix  used  to  fabricate  the  test  block.  The 
accelerometers  and  strain  gage  column  were  then  inserted  into  the  wet 
grout-filled  holes  and  the  displaced  grout  was  removed.  For  the  vertical 
gage  line,  four  strain  gages  were  placed  in  the  second  and  four  in  the 
third  pours. 

The  exterior  concrete  slab  was  carefully  constructed  to  minimize 
shock  attenuation  or  reflection  for  the  air  blast  gages.  After  pouring, 
the  surface  was  polished  to  minimize  surface  roughness.  The  pressure 
gage  canisters  were  installed  with  their  tops  flush  with  the  slab 
surface.  Care  was  taken  to  eliminate  any  surface  irregularities  near  the 
concrete -canister  interface. 

3.5  INSTRUMENTATION 

The  selection  of  airblast  gages  was  based  predominately  on  the 
expected  frequency  and  amplitude  of  the  data.  Sufficient  blast  data 
(discussed  earlier)  were  available  to  predict  both  parameters.  However, 
an  additional  constraint  was  imposed.  It  was  expected  that  dilatation  of 
the  storage  chamber  model  would  impose  a  lateral  acceleration  to  the 
blast  slab.  Thus,  the  pressure  sensors  needed  to  have  negligible  lateral 
acceleration  sensitivity.  As  indicated  in  Table  3.4,  the  manufacturer's 
specifications  indicated  that  the  PCB  pressure  gages  met  these 
requirements.  The  pressure  gages  were  mounted  in  waterproof , 
shock- isolation  canisters.  The  sensing  surface  of  each  gage  was  slightly 
recessed,  The  recessed  region  (approximately  0.06  in.)  was  filled  with 
silicone  rubber  to  minimize  light  and  thermal  sensitivity. 
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The  accelerometers  were  selected  based  on  an  extrapolation  of 
spherical  and  cylindrical  decoupling  data  (discussed  earlier) .  The 
manufacturer's  specifications  for  these  gages  are  presented  in  Table  3.5 
The  accelerometers  were  placed  in  high-strength,  watertight  canisters. 
The  lowest  natural  frequency  of  the  canister  was  experimentally 
determined  to  be  32,000  Hz,  which  was  compatible  with  the  frequency 
response  of  the  recording  system.  The  mounting  base  was  sufficiently 
strong  to  eliminate  any  base-induced  strain  which  wo-’M  degrade  the 
accelerometer  signal. 

To  measure  strain,  eight  waterproofed  strain  gages  were  bonded  to  a 
1.875-inch  diameter  grout  column.  The  strain  gages  were  1/4-inch, 
350-OHM  gages  with  a  gage  factor  of  2.14  at  75  degrees  F.  The  grout 
column  was  made  from  the  same  mix  used  in  the  model  to  insure  matching 
properties.  The  properties  of  the  grout  are  given  in  Table  3.6. 

3 . 6  DESCRIPTION  OF  THE  CHARGE 

The  energy  source  used  in  all  tests  was  Pentereythuital  tetranitrate 
(PETN) ,  with  the  trade  name  Primacord.  PETN  has  several  advantages.  It 
was  readily  available  in  strands  of  various  linear  densities.  As  a 
result,  a  linear  charge  could  be  assembled  easily  by  varying  the  number 
and  the  density  of  the  strands.  The  burn' rate  (26,200  ft/sec)  was  fast 
enough  to  effectively  produce  a  cylindrical  ground  shock.  Additionally, 
experimental  data  and  the  theoretical  techniques  were  available  to 
accurately  predict  the  equilibrium  pressure  developed  in  the  storage 
chamber  by  detonation  of  the  PETN  (Figure  2.2). 

The  charges  were  assembled  by  taping  the  required  number  of  strands 
together.  The  charges  were  detonated  at  the  portal  ends.  The  smallest 
and  largest  charges  used  in  the  test  program  are  shown  in  Figure  3.5. 

3.7  PROCEDURE 

The  following  procedure  was  used  for  each  test; 

a.  The  sensors  were  calibrated. 

b.  The  charge  was  assembled  and  placed  in  the  test  fixture. 

c.  The  tape  machines  were  started. 
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d.  The  charge  was  fired  and  the  data  recorded. 

e.  The  test  fixtures  were  cleaned  for  the  next  test. 

3.8  TEST  PROGRAM 

The  test  program  was  performed  in  three  phases.  Phase  1  was  designed 
primarily  to  obtain  ground  shock  data;  however,  one  or  more  airblast 


stations  were  monitored  during 

following  tests: 

each  test.  Phase  I  consisted  of  the 

i 

Number 

Loading 

Dens i tv 

of  Tests 

lb /ft3 

kg/m3. 

4 

0.1 

1.6 

4 

0.3 

4.8 

4 

1.0 

16.0 

3 

2.0 

32.0 

Several  tests  were  conducted  at  each  loading  density  to  establish  the 
reproducibility  of  the  data.  The  loading  density  is  the  total  charge 
weight  divided  by  the  volume  of  the  storage  chamber.  To  make  sure  that 
the  accelerometers  and  strain  gages  were  appropriately  ranged,  the 
loading  densities  were  gradually  increased  during  the  series. 

Phase  2  was  designed  to  obtain  airblast  pressure  data  outside  the 
tunnel.  It  consisted  of  a  total  of  51  tests  conducted  at  five  loading 
densities,  ranging  between  0.1  and  0.46  lb/ft3.  Low  loading  densities 
were  used  initially  to  prevent  damage  to  the  model.  This  phase  consisted 
of  the  following  tests; 


Number 

Loadine 

Density 

of  Tests 

lb/ft3 

kg/m3 

14 

0.1 

1.6 

13 

0.2 

3.2 

15 

0.28 

4.5 

6 

0.37 

6.0 

3 

0.46 

7.4 

Phase  3  consisted  of  a  single  test  at  a  loading  density  of  25.3 
lb/ft3  (405  kg/m3).  Catastrophic  failure  of  the  model  was  anticipated; 
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however,  peak  strains  were  measured  successfully,  since  the  breakup  and 
mass  motion  of  the  block  occurred  after  passage  of  the  incident  stress 
wave. 

3 . 9  DATA  ACQUISITION  AND  REDUCTION 

The  data  were  recorded  on  a  14- track  FM  analog  tape  machine  operating 
at  a  tape  speed  of  120  in. /sec.  The  frequency  response  of  the  recording 
system  was  40  KHz,  which  was  consistent  with  the  frequency  response  of 
other  system  components.  The  analog  data  were  digitized  on  an 
analog- to-digital  converter,  recorded  on  magnetic  tape,  and  subsequently 
processed  through  the  WES  GE-635  computer.  Velocity  histories  were 
obtained  by  numerical  integration  of  the  acceleration  histories.  Typical 
plots  of  velocity,  strain  and  pressure  histories  are  presented  in 
Figures  3.6,  3.7,  and  3.8,  respectively  (note  that  negative  strains  are 
compressive  on  all  records) .  A  typical  strain  history  obtained  at  a 
loading  density  of  25.3  lb/ft3  is  presented  in  Figure  3.9.  Peak 
velocity,  strain  and  pressure  data  were  obtained  from  these  records  and 
are  recorded  in  Tables  3.7,  3.8,  and  3.9,  respectively.  These  results 
are  analyzed  in  Chapter  4. 
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Table  3.1  Design  Parameters  for  1:75  Scale  Model  and 
Full-Scale  Underground  Storage  Facility 


Parameter 

Model 

Full-Scale 

Access  Tunnel  Length 

3.5  ft 

262.5  ft 

Access  Tunnel  Diameter 

0.223  ft 

16.7  ft 

Storage  Chamber  Diameter 

0.318  ft 

23.8  ft 

Storage  Chamber  Length 

3.76  ft 

282  ft 

Storage  Chamber  Volume 

0.298  ft3 

125,550  ft3 

Charge  Mass  for  Loading 

0.744  lb 

313,875  lb 

Table  3.2  Distance  Between  the  Center  of  the  Storage 
Chamber  and  Ground  Motion  Gages 


Accelerometers 

Strain  Gases 

Distance 

ft 

Gage  Number 

Distance 

ft 

1 

1.65 

8 

1.13 

2 

3.03 

9 

1.46 

4 

4.96 

10 

1.79 

5 

1.38 

11 

2.13 

6 

1.93 

12 

2.46 

7 

(deleted) 

13 

2.79 

14 

3.13 
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Table  3.3  Location  of  Airblast  Measurement  Stations  with  Respect 
to  the  Portal  and  the  Extended  Tunnel  Centerline 


Station  Number 

Angle  with  Respect  to  the 
Extended  Center  line 
degrees 

Radial  Distance 
From  Tunnel  Portal 
ft 

1 

0 

0.75 

2 

0 

1.50 

3 

0 

2.50 

4 

0 

10.00 

5 

30 

3.00 

6 

30 

5.00 

7 

30 

8.00 

8 

45 

2.00 

9 

45 

3.50 

10 

45 

5.00 

11 

45 

7.00 

12 

60 

1.00 

13 

60 

2.50 

14 

60 

5.00 
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Table  3.4  Manufacturer's  Specifications  for  Pressure  Sensors 


Sensitivity 

Resolution 

Resonant  Frequency 

Rise  Time 

Time  Constant 

Low  Frequency  (5%  down) 

Linearity 

Range  (for  5  volts  out) 
Range  (for  10  volts  out) 
Maximum  Pressure 
Output  Impedance 
Vibration/Shock 
Acceleration  Sensitivity 
Temperature  Range 
Temperature  Coefficient 
Flash  Temperature 


50  mV/psi 
0.004  psi 
250,000  Hz 
2  micro- seconds 
1  sec 
1  Hz 

1  percent 
100  psi 
200  psi 
1000  psi 
100  ohms 
2,000/20,000  g 

0.002  psi/g 
-100  to  +270  F 
0.03%  F 
3000  F 


36 


Table  3.5  Manufacturer's  Specifications  for  Accelerometers 


Range,  g  pk 

±5,000 

±10,000 

±20,000 

±50,000 

Sensitivity,  mV/g,  at 

10  Vdc,  Nominal  value: 

Minimum  value: 

0.100 

0.075 

0.050 

0.037 

0.025 

0.018 

0.010 

0.007 

Mounted  Resonance 

Frequency,  Hz,  nominal 

50,000 

70,000 

100,000 

180,000 

Useful  Frequency  Response, 

Hz, 

8,000 

9,000 

12,000 

30,000 

Environmental  Acceleration 
Limits  (Sensitive  Axis,  g  pk) 

±12,500 

±25,000 

±50,000 

±100,000 

Transverse  Axis,  g  pk 

±12,500 

±25,000 

±30,000 

±50,000 

Minimum  Half- Sine  Pulse 
Duration,  microseconds 

125 

90 

65 

30 

Table  3.6  Experimentally-Determined  Grout  Properties 
(average  values  from  six  specimens) 


Specific  Gravity 

2.19 

Compressional  Wave  Velocity  (ft/sec) 

12,954 

Shear  Wave  Velocity  (ft/sec) 

8,373 

Shear  Wave  Velocity/Compressional  Velocity 

0.647 

Young's  Modulus  (psi) 

4.727 

X 

10® 

Shear  Modulus  (psi) 

2.071 

X 

10® 

Bulk  Modulus  (psi) 

2.201 

X 

10® 

Lame  Constant  (psi) 

8.205 

X 

105 

Poisson's  Ratio 


0.140 


Table  3.7  Peak  Velocities  Obtained  as  a  Function  of 
Loading  Density  and  Sange 


Loading 

Density _ Peak  Particle  Velocity,  ft/sec 


Table  3.8  Peak  Strain  Obtained  as  a  Function  of  Loading  Density 
and  Range 


Loading 

Densitv 

Peak  Strain, 

,  yin. /in. 

lb/ft3 

"$'792 

1,125 

1,458 

Range  (ft) ' 

1,792  ZA15 

2,45.8 

2J92 

3,125 

0.1 

54 

69 

38 

40 

33 

23 

23 

16 

0.1 

40 

78 

25 

— 

.... 

.... 

.... 

0.1 

40 

107 

30 

37 

20 

14 

.... 

.... 

0.1 

34 

20 

20 

12 

18 

13 

8 

8 

0.3 

119 

.... 

86 

49 

52 

53 

22 

12 

0.3 

97 

94 

43 

33 

33 

27 

25 

.... 

0.3 

119 

109 

102 

63 

61 

54 

41 

33 

0.3 

139 

130 

95 

57 

63 

55 

44 

48 

1.0 

373 

224 

261 

■  <•  m 

71 

41 

39 

28 

1.(1) 

240 

229 

155 

»  w  m  m 

54 

42 

24 

26 

1.0 

310 

271 

206 

44 

32 

22 

30 

21 

1.0 

347 

237 

250 

47 

43 

58 

25 

31 

2.0 

560 

626 

465 

224 

149 

86 

70 

68 

2.0 

772 

.... 

633 

291 

290 

145 

137 

63 

2.0 

984 

956 

767 

528 

-  »*•  • 

•»  • 

as.  3 

»  •  «*  • 

7000 

7100 

7800 

4985 

3850 

2420 

25.3 

•  m  m  » 

9000* 

5800* 

5900* 

4800* 

W>  U  M  H 

.... 

*  Obtained  from  a  second  strain  gage  column  used  only  for  that  test. 


Table  3.9a  Peak  Blast  Pressures  Measured  Along  the 
0- Degree  Line  (Extended  Centerline  of 
the  Access  Tunnel) 


ing  Density 
lb /ft3 

Radial  Distance 
ft 

Peak  Pressure 
DSi 

0.10 

0.75 

25.0 

0.10 

0.75 

23.0 

0.10 

0.75 

28.0 

0.10 

0.75 

23.7 

0.10 

10.0 

0.21 

0.10 

10.0 

0.27 

0.10 

10.0 

0.22 

0.10 

10.0 

0.23 

0.19 

10.0 

0.81 

0.19 

10.0 

0.46 

0.19 

10.0 

0.29 

0.19 

10.0 

0.32 

0.28 

10.0 

1.31 

0.28 

10.0 

0.22 

0.28 

10.0 

1.05 

0.28 

0.75 

43.0 

0.28 

0.75 

46.0 

0.28 

0.75 

48.0 

0.28 

0.75 

40.0 

0.97 

0.75 

70.0 

0.97 

1.50 

37.5 

0.97 

10.0 

3.76 

0.97 

1.50 

40.0 

0.97 

2.50 

24.1 

0.97 

2.50 

18.0 

1.95 

10.0 

6,6 

1.95 

10.0 

6.2 

1.95 

2.5 

55.0 

1.95 

10.0 

6.8 
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Table  3.9b  Peak  Blast  Pressures  Measured  Along  the 
30 -Degree  Line  (Measured  from  Extended 
Centerline  of  the  Access  Tunnel) 


Loading  Density 
lb /ft3 

Radial  Distance 
ft 

Peak  Pressure 
DSi 

0.10 

3 

0.91 

0.10 

3 

0.75 

0.10 

5 

0.82 

0.10 

5 

0.65 

0.10 

8 

0.63 

0.10 

8 

0.68 

0.10 

8 

0.30 

0.19 

3 

1.7 

0.19 

3 

1.9 

0.19 

5 

1.4 

0.19 

5 

1,57 

0.19 

8 

0.70 

0.19 

8 

0.69 

0.28 

3 

4.35 

0.28 

3 

4.9 

0.28 

5 

1.5 

0.28 

5 

1.5 

0.28 

8 

0.85 

0,28 

8 

0.78 

0.37 

3 

4.8 

0.37 

3 

5.8 

0.37 

5 

2.38 

0.37 

5 

2.4 

0.37 

.  8 

1.6 

0.37 

8 

1.59 

0.46 

3 

8.4 

0.46 

5 

3.6 

0.46 

8 

3.0 
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Table  3.9c  Peak  Blast  Pressures  Measured  Along  a 
45 -Degree  Line  with  Respect  to  the  Ex¬ 
tended  Center  line  of  the  Access  Tunnel 


Loading  Density 
lb /ft3 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

0.19 

0.19 

0.19 

0.19 

0.19 

0.19 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

0.37 

0.37 

0.37 

0.37 

0.37 

0.37 

0.37 

0.37 

0.46 

0.46 

0.46 

0.46 


Radial  Distance 
ft _ 


2.0 

2.0 

3.5 

3.5 

3.5 

7.0 

7.0 

7.0 

2.0 

2.0 

3.5 

3.5 

7.0 

7.0 

2.0 

2.0 

3.5 

3.5 

5.0 

5.0 

7.0 

7.0 

2.0 

2.0 

3.5 

3.5 

5.0 

5.0 

7.0 

7.0 

2.0 

3.5 

5.0 

7.0 


Peak  Pressure 
nsi 

1.0 

0.9 

1.0 

0.75 

0.64 

0.58 

0.40 

0.40 

1.88 

2.25 
1.9 
2.08 
0.71 
0.69 

8.3 

7.5 
3.0 

2.6 
1.07 
1.08 
0.8 
0.75 

11.2 

7.3 
3.98 

3.85 

1.86 
1.6 
1.56 

1.25 

12.5 

5.3 

3.3 
3.1 
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Table  3.9d  Peak  Blast  Pressures  Measured  Along  the 

60 -Degree  Line  (Measured  from  the  Extended 
Centerline  of  the  Access  Tunnel) 


Loading  Density 
lb /ft3 

Radial  Distance 
ft 

Peak  Pressure 
DSi 

0.10 

1.0 

5.0 

0.10 

.  l.o 

4.25 

0.10 

2.5 

0.65 

0.10 

2.5 

0.55 

0.10 

5.0 

0.70 

0.10 

5.0 

0.68 

0.10 

5.0 

0.68 

0.19 

1.0 

9.8 

0.19 

1.0 

12.5 

0.19 

2.5 

1.08 

0.19 

2.5 

1.1 

0.19 

5.0 

0.92 

0.19 

5.0 

0.80 

0.28 

1.0 

34.5 

0.28 

1.0 

36.5 

0.28 

2.5 

3.85 

0.28 

2.5 

3.9 

0.28 

5.0 

1.0 

0.28 

5.0 

l.o 

0.37 

2.5 

4.68 

0.37 

2.5 

3.7 

0.37 

5.0 

1.5 

0.37 

5.0 

1.54 

0.46 

2.5 

5.25 

0.46 

5.0 

2.6 
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Figure  3.2.  Schematic  indicating  the  location  of  the  ground  shock  gages 


Figure  3. 


.  Layout  (plan  view)  of  external  airblast 
gages  referenced  to  tunnel  portal. 


47 


(a)  End  view  after  first  pour  (b)  Side  view  after  first  pour 


(d)  Completed  model 


Figure  3.4,  Construction  of  the  Small  Scale  Model 


Figure  3.5.  The  Smallest  and  Largest  Charge  Used  in  the  Test  Program. 


ACCELERATION 


i  i  i - 1 - r - 1 

0.000  0.001  0.002  0.003  0.004  0.005 

TIME  FROM  GET  -  SEC 


Figure  3.6.  Acceleration  history  obtained  at  a  loading  density  of 
1  lb/ft^.  Gage  range  is  1.378  ft  from  chamber  center 
line. 
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MICRO  IN  PER  IN 


0.000  0.002  0.004  0.006  0.008  0.010 


TIME  FROH  OET  -  SEC 


Figure  3.7.  Free-fiold  strain  measured  at  a  loading  density  of  1  lb/ft^. 
Cage  range  is  1.458  ft  from  chamber  center  line. 
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TIME  IfiOH  DEI  -  SEC 


figure  3.8.  Free-field  pressure  history  obtained  at  a  loading 
density  of  2  lb/i'tA  Gage  on  0-degree  line 
2.5  ft  froo  tunnel  portal. 


CHAPTER  4 


ANALYSIS  OF  DATA 

4.1  INTRODUCTION 

The  ground  shock  and  airblast  data  were  analyzed  and  compared  to 
previously  published  (through  1978)  data.  Normalized  plots  were 
developed  for  the  ground  shock  data  which  identify  the  point  of 
transition  from  cylindrical  to  spherical  decay,  Dimensionless  equations 
were  developed  which  define  the  upper  bound  of  the  peak  free- field 
strains  produced  by  an  explosion  in  a  cylindrical  storage  chamber  of 
arbitrary  size. 

4.2  GROUND  SHOCK 

The  peak  particle  velocities  (v)  obtained  in  this  investigation  are 
plotted  in  Figure  4.1  as  a  function  of  r/a,  where  "r"  is  the  radial 
distance  from  the  center  of  the  cavity  to  the  point  of  measurement  and 
"a"  is  the  radius  of  the  storage  chamber.  The  scatter  is  consistent  with 
that  obtained  in  typical  ground  shock  investigations.  The  peak  velocity 
data  agrees  well  with  data  (Amend,  References  21  and  22)  obtained  from 
explosions  in  2-ft  diameter  cavities  in  sandstone  and  highly- jointed 
granite.  These  data  represent  an  order  of  magnitude  increase  in  cavity 
diameter,  however. 

The  peak  velocity  is  strongly  dependent  upon  loading  density. 
Comparing  the  data  obtained  at  0.1  and  1  lb/ft3,  one  finds  that  an  order 
of  magnitude  increase  in  loading  density  results  in  approximately  an 
order  of  magnitude  increase  in  peak  velocity. 

Peak  strain  data  were  derived  from  the  peak  velocity  data  by  using 
the  approximation 

e  ••  v/c  (4.1) 

where  «  is  peak  strain 

v  is  peak  particle  velocity  (in. /sec) 
and  c  is  compressions!  wave  velocity  (in. /sec) 

The  peak  strains  calculated  from  the  peak  velocity  measurements  are 
compared  with  the  measured  peak  strains  (Table  3.7)  In  Figures  4.2 
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through  4.6.  To  facilitate  interpretation,  the  raw  data  obtained  for 
each  loading  density  are  presented  in  separate  figures.  The  scatter  in 
the  data  is  consistent  with  the  scatter  typically  obtained  from  strain  or 
velocity  measurements.  It  should  be  noted  that  the  strain  measurements 
were  obtained  from  a  vertical  array  of  gages,  whereas  the  velocity  data 
was  obtained  from  a  horizontal  array.  Agreement  of  the  data  demonstrates 
that  there  were  no  apparent  directional  anomalies. 

The  data  obtained  near  the  chamber  is  consistent  with  previously 
published  results  obtained  at  radial  distances  equal  to  several  cavity 
lengths  (Atchison,  Reference  19).  This  is  indicated  in  Figure  4.6  by 
comparing  the  data  at  the  25.3  lb/ft3  loading  with  Atchison's  data 
obtained  at  21.6  lb/ft3. 

As  demonstrated  by  Drake  (Reference  20),  Atchison's  data  is  dominated 
by  spherical  decay.  As  indicated  in  Figure  4.6,  the  slope  of  the  data 
obtained  in  this  investigation  appears  to  asymptotically  approach  a 
spherical  decay,  although  the  precise  point  of  transition  from 
cylindrical  to  spherical  decay  cannot  be  clearly  ascertained.  It  can, 
however,  be  determined  from  a  normalized  plot  as  a  result  of  an  increase 
in  data  density. 

The  strain  at  the  cavity  wall  is  determined  essentially  by  the 
chamber  pressure  and  the  properties  of  the  medium.  As  a  first 
approximation,  the  strain  at  the  wall  is  given  by 

ew  "  V*o  (4*2) 

where  ew  is  the  strain  at  the  wall 

P0  is  the  chamber  pressure  (psi) 
and  E0  is  the  elastic  modulus  of  the  medium  (psi) 

The  chamber  pressure  is  determined  by  the  loading  density  and  the 
chemical  composition  of  the  explosive.  It  may  be  determined  for  PETN 
from  Figure  2.2,  Since  the  strain  at  the  wall  is  the  maximum  strain, 
Equation  4.2  may  be  employed  to  normalize  the  preceding  strain  data. 

The  normalized  data  presented  in  Figure  4.7  contains  strain  data  for 
loading  densities  encompassing  over  two  orders  of  magnitude  (0.1  to 
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25.3  lb/ft3).  The  data  coalesces  with  no  increase  in  data  scatter. 
Excellent  agreement  exists  between  the  data  obtained  in  this 
investigation  and  both  the  near-field  data  reported  by  Amend 
(References  21  and  22)  and  the  far -field  data  reported  by  Atchison, 
et  al.  (Reference  19). 

The  transition  from  cylindrical  to  spherical  decay  occurs  at  a  value 
of  r/a  approximately  equal  to  11.8,  which  corresponds  to  a  radial 
distance  equal  to  one -half  the  length  of  the  cavity.  The  rate  of  decay 
of  peak  strain  in  the  cylindrical  region  where  -  <  ,  and  in  the 

spherical  region,  where  is  r"1/2  and  r'2,  respectively.  Although 

the  cylindrical  decay  region  is  apparent,  data  were  not  obtained  in  this 
investigation  at  radial  distances  less  than  L/4  (or  r/a  <  5) .  The  data 
reported  by  Amend  for  |  -  3  are  included.  The  line  which  appears  to 
best  fit  the  data  in  the  cylindrical  decay  region  is  extrapolated  to  the 
cavity  wall  (r/a  *»  1)  in  Figure  4.7.  The  theoretical  intercept  should  be 
1.0  compared  to  the  suggested  value  of  1.2.  Considering  the  range  of 
loading  densities  investigated,  inelastic  effects  near  the  cavity 
surface,  and  the  scatter  inherent  in  strain  data,  the  experimentally 
implied  intercept  agrees  remarkably  well  with  theory. 

Lines  which  identify  the  upper  bound  of  the  data  obtained  from  this 
particular  geometry  are  indicated  in  Figure  4.7.  Employing  the  above 
results,  the  upper  bound  of  the  peak  free- field  strain  produced  by  an 
explosion  in  an  arbitrary  cylindrical  cavity  can  be  determined.  It  has 
been  established  in  this  investigation  that  the  transition  from 
cylindrical  to  spherical  decay  will  occur  at  radial  distances  on  the 
order  of  one -half  the  cavity  length.  The  slopes  in  the  respective  region 
have  been  demonstrated  in  this  and  several  other  investigations. 
Additionally,  the  upper  bound  of  the  normalized  strain  intercept  was 
determined  to  be  2.4,  Using  these  results,  the  upper  bound  equations 
were  found  to  be 

eE0/po  ■  2.4(r/a)“l/2.  when  r/a  <  L/2a  (4.3a) 

and 


eE0/P0  ■  0,85(L/a) l*5(r/a)”2;  when  r/a  >  L/2a  (4.3b) 
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For  larger  standoff  distances,  a  more  convenient  equation  is 
desirable.  From  Figure  2.1,  it  can  be  observed  that  the  chamber  pressure 
generated  by  a  TNT  explosion  varies  nearly  linearly  with  loading  density 
for  loading  densities  between  0.1  and  25.3  lb/ft3.  As  a  result,  the 
chamber  pressure  can  be  closely  approximated  by 

P0  -  25  80  (  7)0,981 

where  P0  is  the  chamber  pressure  (psi)  and  7  is  the  TNT -equivalent 
loading  density  (lb/ft3).  However,  to  simplify  the  equation  which  will 
be  developed  and  to  increase  conservatism,  an  exponent  of  unity  is 
assumed . 

Substituting  for  P0  and  rearranging  Equation  4.3  we  obtain 
r/a  =  3.83  x  10^(y/eEo)2;  when  r/a  L/2a  (4.4a) 

and 

r/a  =  46.8(L/a)0‘75(Y/eEo)l/2;  when  r/a  L/2a  (4.4b) 

Equation  4.4  is  employed  in  the  next  chapter  to  predict  the  safe  standoff 
distances  required  for  ground  shock.  Since  the  equation  predicts  the 
maximum  distance  at  which  a  specified  strain  will  occur,  the  predicted 
standoff  distances  will  be  safely  conservative. 

4.3  AIRBLAST 

The  peak  free -field  airblast  data  obtained  in  this  investigation  are  . 
tabulated  in  Table  3.8.  These  data  are  presented  in  Figures  4.8  through 
4.11  as  a  function  of  the  exit  pressure  (PE),  diameter  of  the  access 
tunnel  (D) ,  and  range  (r) .  The  exit  pressure  at  the  portal  was 
determined  from  Equation  2.1.  The  nondimensionalized  peak  free-field 
blast  data  obtained  along  the  extended  centerline  and  the  30,  45  and 
60-degree  lines  (with  the  0-degree  line  being  the  extended  centerline) 
are  presented  in  Figures  4.8  through  4.11,  respectively.  Data  obtained 
by  other  investigators  are  also  included  and  upper  bound  curves  are 
indicated. 
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Data  obtained  along  the  extended  centerline  (Figure  4.2)  with  loading 
densities  of  0.97  and  1.95  lb/ft3  agree  with  data  obtained  by  other 
investigators  (Skjeltorp,  References  4  and  10)  for  exit  geometries  with 
0  and  30-degree  surface  slopes.  This  result  indicates  that  the  standoff 
distance  required  along  the  extended  centerline  is  not  sensitive  to  exit 
slopes  between  0  and  30  degrees.  At  low  loading  densities  (0.1,  0.19  and 
0.28  lb/ft3),  the  data  obtained  at  large  scaled  distances  (r/D  -  44.8) 
had  substantial  scatter.  This  scatter  does  not  appear  to  be 
geometrically  dependent,  however. 

Substantial  data  scatter  exists  in  all  of  the  data  obtained  along  the 
30,  45  and  60-degree  lines  (Figures  4.9,  4.10  and  4.11).  It  is  the 
author's  opinion  that  this  occurred  because  these  experiments  were 
conducted  with  low  loading  densities  (0.1  to  0.46  lb/ft3)  to  prevent 
damage  to  the  model  prior  to  completion  of  the  ground  shock 
investigations.  As  discussed  earlier,  low- strength  shocks  produced  by 
low  loading  densities  appear  to  be  substantially  attenuated  by 
microscopic  wall  roughness.  This  results  in  lower  exit  pressures  than 
predicted  and,  consequently,  a  reduced  level  of  shock  propagates  outside 
the  portal.  These  effects  are  nonlinear ly  related  to  the  loading 
density.  Since  the  loading  density  is  used  to  calculate  the  exit 
pressure,  data  scatter  results.  Gurke  and  Scheklinski  (Reference  15) 
reported  similar  results  at  low  loading  densities. 

Based  upon  the  data  presented  in  Figures  4.8  through  4.11,  the  upper 
bound  of  the  free-fiald  blast  pressure  can  be  expressed  in  the 
dimensionless  form: 

P/PE  =  Cq(R/D)”3/2  (4.5) 

where  R  is  the  radial  distance  from  the  portal  (ft) 

D  is  the  effective  diameter  of  the  access  tunnel  (ft) 

P  is  the  peak  blast  pressure  (psi)  at  R 
PE  is  the  peak  overpressure  at  the  tunnel  exit  (psi) 
and  Cq  is  a  directional  coefficient  dependent  upon  the  direction  of 
measurement  with  respect  to  the  extended  center  line. 
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The  values  of  Cq  as  a  function  of  0  are: 

o  o  o  o 

0:  0  30  45  60 

C0:  3.8  3.2  3.2  2.0 

Equation  4.5  is  consistent  with  Equation  2.3  developed  by  Skjeltorp 
(Reference  4).  The  dimensionless  form  of  Equation  4.5  provides  a  method 
for  predicting  the  standoff  distance  required  as  a  function  of  a  blast 
pressure  threshold  for  structures  in  the  particular  area  of  interest,  or 
as  a  function  of  different  peak  pressure  criteria  used  in  various 
countries. 

In  all  directions  from  the  portal,  the  upper  bound  of  the  airblast 
data  obtained  in  this  investigation  is  consistent  with  the  results 
obtained  by  other  investigators.  Also,  the  rate  of  peak  pressure  decay 
is  consistent  with  Westine's  (Reference  16)  gun  blast  data  and  the 
free-field  blast  data  obtained  by  other  investigators. 

4.4  SUMMARY 

Equations  have  been  developed  which  predict  the  upper  bound  of  the 
free-field  strain  produced  by  decoupled  explosions  in  cylindrical 
cavities  of  finite  length,  sited  in  competent  rock.  Cylindrical  wave 
propagation  was  found  to  be  dominant  out  to  a  maximum  radial  distance 
(from  the  center  of  the  cavity)  equal  to  one -half  of  the  cavity  length. 
At  larger  radial  distances,  the  peak  free-field  parameters  will  be 
dominated  by  spherical  wave  propagation.  These  results  may  be  used  to 
determine  safe  standoff  distances  related  to  ground  shock  effects.  The 
upper  bound  of  the  airblast  data  obtained  in  this  investigation  was 
consistent  with  that  found  by  other  investigators. 
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Figure  4.1.  Peak  particle  velocity  versus  dimensionless  distance 
(distance  "R"  divided  by  chamber  radius  WAM)  from 
the  center  of  cavity. 
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Figure  4.2.  Peak  strain  as  a  function  of  dimensionless  distance 

(distance  "ft”  divided  by  chamber  radius  "A")  obtained 
at  a  loading  deusity  of  0.1  lb/ftA 
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Figure  4.3.  Peak  strain  as  a  function  of  dioenaionlees  distance 

(distance  “R1*  divided  by  chamber  radius  '‘A'1)  obtained 
at  a  loading  density  of  0.3  lb/ft^. 
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Figure  4.4. 


Peak  strain  as  a  function  of  distension  less  distance 
(distance  '’R'*  divided  by  cbaober  radius  "AM)  obtained 
at  a  loading  density  of  1.0  lb/ft^. 
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Figure  4.5.  Peak  strain  as  a  function  of  dimensionless  distance 

(distance  "R"  divided  by  chamber  radius  MA")  obtained 
at  a  loading  density  of  2.0  lb/ft^. 
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Figure  4.8.  Normalized  blast  data  (measured  peak  pressure  MP°  divided 
by  chamber  pressure  "Pe")  along  the  extended  centerline 
as  a  function  of  dimensionless  distance  (distance  **8“ 

divided  by  the  access  tunnel  diameter 
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Figure  4.9.  Normalized  blast  data  (measured  peak  pressure  MPM  divided 
by  chamber  pressure  "Pe")  along  a  30-degree  line  as  a 
function  of  dimensionless  distance  (distance  MR"  divided 
by  access  tunnel  diameter  "D") . 
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All  data  from 
WES  tests 


bigute  4.10.  Normalized  blast  data  (measured  peak  pressure  "PH  divided 
by  chamber  pressure  MPe”)  along  a  45-degree  line  as  a 
function  of  dimensionless  distance  (distance  "R"  divided 
by  access  tunnel  diameter  "D"). 
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Figure  4.11.  Normalised  blast  data  (measured  peak  pressure  "P" 
divided  by  chamber  pressure  along  a  60-degree 

line  as  a  function  of  dimensionless  distance 
(distance  “S'*  divided  by  access  tunnel  diameter  "D") . 
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CHAPTER  5 


PROPOSED  EQUATIONS  FOR 
SAFETY  STANDARDS 


5.1  INTRODUCTION 

Based  on  the  results  of  this  study,  equations  were  developed  to 
determine  the  minimum  safe  spacing  of  underground  storage  chambers  to 
prevent  explosive  communication  or  damage  to  stores  in  adjacent  chambers. 
The  analysis  assumes  that  the  chambers  are  sited  in  competent  rock. 

These  equations  are  based  on  the  peak  free -field  strain  produced  by 
decoupled  explosions  in  cylindrical  cavities.  The  standoff  distance 
(from  the  storage  chamber)  required  to  prevent  ground  shock  damage  to 
inhabited  buildings  is  a  function  of  the  maximum  particle  velocity  that 
may  be  safely  induced  at  the  base  of  the  structure  at  risk.  Predictive 
techniques  to  determine  the  safe  airblast  standoff  distance  are  also 
reviewed.  The  safety  equations  determined  in  this  investigation  are 
compared  to  those  currently  employed  in  Reference  1,  "DOD  Ammunition  and 
Explosive  Safety  Standards*1  (AESS),  .  v  • 

5 . 2  MINIMUM  CHAMBER  SEPARATION  REQUIRED  TO  PREVENT,  EXPLOSIVE  PROPAGATION 

The  distance  between  underground  storage  chambers  must  be  large  /  / 

enough  to  prevent  oxplosive.  communication  as  a  result  of  catastrophic  * 
wall  failure  or  spalling.  It  Was  demonstrated  in  Chapter  2  that  ' 

catastrophic  wall  failure  will  occur  in  unlinodtuntxels  if  the  peak  ; ' 
strain  at,  the  wall  is  greater  than  0.004,  or:  0.4"pereent.  Equation  4.4, 
which  defines  the  relationship  between  radial  distance  r,  peak  free -field 
strain  e,  TNT- equivalent  lojidingdonsity  y,  the  Radius  of  the  cavity  a, 
and  the  properties  of  the  medium^ ^can  he  used; to  determine  the  required 
spacing.  To  apply  this  equation  to  noncyi indr teal  cavities,  an  effective 
radius  must  be  employed.  The  effective  radius  is  defined  by 

ati  “  <A/4n)  w  •  (5-1) 

whore  a.  is  the  effective  radius  in  ft 

and  A  is  the  cross-sectional  area  of  the  storage  chamber  in  ft2. 
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Substituting  the  strain  threshold  for  catastrophic  failure  into 
Equation  4.4,  we  obtain 

rc/ae  =  2.39  x  1012(y/Eo)2;  when  rc/ae  <  L/2ae  (5.2a) 

and 


rc/ae  =  7.4  x  102(L/ae)^* 75(y/EQ) 1/2.  yhen  rc/ae  >L/2ae  (5.2b) 

where  r0  is  the  standoff  distance  (in  feet)  required  to  prevent  tunnel 
closure . 

Similarly,  the  chamber  spacing  required  to  prevent  spall -induced, 
sympathetic  detonations  can  be  obtained  from  Equation  4.4,  if  the  spall 
velocity,  v,p,  required  to  initiate  the  contents  of  the  adjacent  chamber 
is  known  and  can  be  expressed  as  a  function  of  the  peak  free-field 
strain.  For  one - dimens ional  wave  propagation,  particle  velocity  vp  is 
related  to  free-field  strain  by 

vp  -  cp  c  (5.3) 

where  c  is  the  compressional  wave  velocity. 

Although  not  strictly  valid,  Equation  5.3  can  be  used  (as  a  first 
approximation)  for  both  cylindrical  and  spherical  wave  propagation.  As 
indicated  earlier,  the  maximum  spall  velocity  that  can  be  generated  at  a 
free  surface  is  equal  to  twice  the  incident  free-field  particle  velocity. 
Nonce,  the  relationship  between  strain  and  maximum  spall  velocity  is 

tSp  -  vsp/2c  (5.4) 

where  v#p  is  the  upper  bound  of  spall  velocity  at  the  wall. 

Substituting  for  strain  into  Equation  4,4,  the  chamber  separation 
(r»p)  required  to  prevent  spall  from  initiating  explosions  In  adjacent 
magazines  is  given  by 

rap/ae  -  1.53  x  iO^cy/v^Eo)2;  when  rsp/ae<L/2a0  (5.5a) 
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.  (5.5b) 


rsp/ae  =  66.2(L/ae)°*75(cY/vspE0)°‘5;  when  rsp/ae  >  L/2ae 


At  low  loading  densities,  the  chamber  spacing  will  be  determined  by  5.2a 
or  5.5a,  which  apply  to  the  cylindrical  decay  region.  As  the  loading 
density  is  increased,  the  required  spacing  will  increase.  Above  some 
threshold  loading  density,  the  chamber  spacing  is  determined  by 
Equation  5 . 2b  and  5 . 5b  (which  are  based  upon  spherical  decay) .  In  this 
region,  the  standoff  distance  increases  with  an  increase  in  loading 
density,  chamber  length,  or  chamber  radius. 

The  chamber  spacing  required  to  prevent  explosion  communication  is 
determined  by  the  larger  of  the  two  values  (rc  or  rBp) .  However,  it  is 
not  necessary  to  evaluate  both  r0  and  rap  to  determine  the  required 
chamber  spacing.  If  «p  >  0.004,  catastrophic  failure  (Equation  5.2) 
dictates  the  spacing;  conversely  if  ep  <  0.004,  spalling  (Equation  5.5) 
determines  the  chamber  spacing. 

In  most  cases,  the  minimum  chamber  spacing  will  be  determined  by 
catastrophic  wall  failure.  This  is  because  the  spall  velocity  required 
to  initiate  explosives  is  typically  on  the  order  of  400  ft/sec,  which 
corresponds  to  a  strain  at  the  wall  that  is  well  above  the  threshold 
strain  (0.004)  required  for  catastrophic  failure.  Spalling  dominates 
only  when  sensitive  explosives  (initiated  by  very  low  spall  velocities) 
are  in  the  adjacent  chamber.  For  example,  for  granite  (assuming  c  equal 
to  15,000  ft/soc),  spalling  dominates  only  if  the  impact  sensitivity  of 
the  contents  is  loss  than  120  ft/sec. 

The  equation  used  in  the  AESS  (Reference  1)  to  determine  the  chamber 
separation  distance  (rcoo)  required  to  prevent  explosive  communication  is 

tea,  -  1.5  Wi/a  (5.6) 


where  V  is  the  total  explosive  weight  in  the  magazine  (lbs). 

It  is  difficult  to  make  a  general  comparison  between  the  equations 
developed  above  and  Equation  5.6.  because  the  latter  does  not  explicitly 
address  geometrical  effects,  spalling  velocity,  properties  of  the  medium. 
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and  nonlinear  effects  of  decoupling.  However,  the  equations  can  be 
compared  by  assuming  some  geometrical  design,  loading  density,  and  spall 
velocity.  For  comparison  purposes,  the  following  values  were  assumed: 


Chamber  length: 

Chamber  radius: 

Spall  velocity  required  for 
initiation  of  contents: 

Wave  speed  of  rock: 

Young's  Modulus: 


262  ft 
12.24  ft 
400  ft/sec 


15,000  ft/sec 
5  x  106  psi 


These  dimensions  are  representative  of  a  full* scale  facility  design.  The 
distances  required  between  the  chambers,  as  predicted  by  the  above 
equations,  for  spalling  (rip),  tunnel  closure  (r0),  and  the  current  safety 
standards  (rcoa)  are  presented  in  Figure  5.1  for  loading  densities  between 
0.1  and  25  lb/ft3. 

The  values  shown  in  Figure  5.1  indicate  that  the  current  safety 
standards  for  chamber  separation  distance  required  to  prevent  explosive 
communication  are  very  conservative.  As  indicated  earlier,  spalling 
should  normally  present  no  problem.  In  this  example,  a  loading  density 
of  more  than  11  lb/ft3  is  required  to  generate  a  spall  velocity  of 
400  ft/sec  at  the  tunnel  wall.  According  to  Equation  5.2,  zero  chamber 
separation  is  required  to  prevent  tunnel  closure  if  the  loading  density 
is  below  3  lb/ft3.  This  is  because  explosions  in  cavities  with  loading 
densities  less  than  3  lb/ft3  produce  strains  at  the  cavity  wall  that  are 
less  than  the  incident  strain  (0.004)  normally  required  to  produce  wall 
failure.  Based  on  structural  considerations  alone,  a  chamber  separation 
of  one  cavity  diameter  is  recommended.  Loading  densities  on  the  order  of 
2  to  3  lb/ft3  are  currently  being  considered.  At  these  loading 
densities,  standard  construction  procedures  will  preclude  explosive 
communication  in  competent  rock. 

The  above  example  demonstrates  that  current  safety  standards  may 
significantly  overestimate  the  chamber  spacing  required  to  prevent 
explosive  communication.  This  occurs  because  the  current  standards 
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involve  a  lumped  approach,  and  do  not  take  into  consideration  geometrical 
effects,  cylindrical  geometry,  decoupling,  or  the  properties  of  the 
medium. 

5.3  MINIMUM  CHAMBER  SEPARATION  REQUIRED  TO  PREVENT  DAMAGE  TO  STORES 

An  underground  explosion  will  propagate  a  strain  pulse  through  the 
surrounding  rock,  toward  the  wall  of  an  adjacent  (acceptor)  storage 
chamber.  Damage  to  the  ammunition  contents  may  occur  if  the  strain  at 
the  acceptor  chamber  wall  is  sufficient  to  induce  intermittent  damage  to 
the  wall.  Based  upon  experimental  data  obtained  from  several  different 
media,  intermittent  boundary  failures  will  not  occur  at  peak  strains  less 
than  0.0003.  Substituting  for  strain  in  Equation  4.4,  the  minimum 
spacing  (rcd)  required  between  storage  chambers  to  prevent  damage  to  the 
contents  is  given  by 

rcd^ae  "  2,39  x  10U(y/Eo>2  J  when  *c<j/ae  <  L/2ae  (5.7a) 

rcd/ae  ■  2.34  x  103(L/a)0,75(Y/Eo)1^2  ;  when  *C(j/ae  >  h/2aQ  (5.7b) 

According  to  the  AESS ,  the  chamber  separation  (Ded)  required  to 
prevent  damage  to  stored  ammunition  can  be  calculated  by  the  following 
formulas; 

Dgd  -  3.5  W1/3  (sandstone)  (5.8a) 

Dcd  -  4.3  «w  (limestone)  (5.8b) 

Dcd  -  5.0  U1'3  (granite)  (5.8c) 

whore  W  is  the  weight  of  the  explosive  contents  in  lbs. 

As  in  the  previous  case,  a  general  comparison  of  Equations  5,7  and 
5.8  Is  difficult;  however,  they  were  compared  here  using  the  geometrical 
design  and  medium  properties  used  in  the  previous  example.  The  chamber 
spacings  required  by  AESS  (Dsd)  and  by  Equation  5.7  (rcd)  are  presented  in 
Figure  5.2  for  loading  densities  from  0*1  to  25  lb/£t. 

Theoretically,  zero  chamber  separation  is  required  for  loading 
densities  below  0.3  lb/ft.  Explosions  in  storage  chambers  with  low 


loading  densities  result  in  wall  strains  below  the  incident  peak  strain 
typically  required  to  damage  unlined  tunnels  in  rock.  For  the  same 
reasons  discussed  in  the  preceding  example,  a  minimum  chamber  separation 
of  one  cavity  diameter  is  recommended. 

For  this  particular  geometrical  configuration,  the  current  safety 
standards  are  conservative.  At  loading  densities  of  2  and  3  lb/ft3, 
Equation  5.7  suggests  a  chamber  separation  of  167  and  206  ft, 
respectively,  as  compared  to  the  AESS  requirement  of  328  and  375  ft, 
respectively.  The  AESS  standards  are  conservative  for  most  magazine 
geometries.  Equation  5.7  should  more  accurately  predict  the  required 
chamber  spacing. 

5.4  INHABITED  BUILDING  STANDOFF  DISTANCE 

The  required  separation  between  the  storage  chamber  and  inhabited 
buildings  can  be  determined  from  Equation  4.4  by  employing  the 
approximation  given  in  Equation  5.3.  Substituting  for  strain  into 
Equation  4.4: 

rib/ae  ■  x  107(cYyvibEo)2i  when  rib/ae  L/2ae  (5.9a) 


and 

*ib/ac  "  8  (L/ae)8*73(cY/VibE0)8,3;  when  rj.b/ae  L/2ae  (5,9b) 

where  rlb  is  the  standoff  distance  (ft)  required  for  inhabited  buildings, 
and  Vtb  is  the  maximum  ground  motion  velocity  (ft/sec)  which  may  be 
safely  induced  at  the  base  of  the  structures. 

Equation  5.9  defines  the  distance  from  the  center  of  a  storage 
chamber  required  for  the  particle  velocities  to  decay  to  a  specified 
value  (Vu,),  as  a  function  of  the  loading  density,  dimensions  of  the 
cavity,  and  properties  of  the  medium.  The  equations  used  in  AESS  do  not 
address  many  of  these  parameters.  According  to  AESS,  the  standoff 
distance  required  for  inhabited  buildings  can  be  determined  by: 

Djj,  »  2, 33(Y/tt)8* 3  «4/9  (saaii,  gravel,  moist  clay)  (5.10a) 

0^  -  U.8(Y/Yt)0*3  W^/9  (soft  rock)  (5.10b) 
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Dib  -  13.3(7/7t)0-3  W*/9  (hard  rock) 


(5.10c) 


where  Dlb  Is  the  standoff  distance  (ft) 

W  is  the  total  TNT -equivalent  charge  weight  (lbs) 

7  is  the  loading  density  (lb/ft3) 

7fc  is  the  density  of  TNT  (lb/ft3) 

These  equations  are  predicated  upon  maximum  allowable  particle  velocities 
of  2.4,  4.5  and  9  in. /sec,  respectively,  at  the  base  of  structures  in 
these  geologies. 

The  AESS  equations  have  several  disadvantages.  They  are  apparently 
based  upon  the  decoupling  effects  of  spherical  cavities  and  consequently 
do  not  take  into  consideration  the  additional  decoupling' effects  of 
cylindrical  cavities  of  finite  length.  The  predictive  technique  also 
does  not  offer  the  capability  to  calculate  the  standoff  distance  required 
tor  structures  that  might  be  damaged  at  significantly  higher  or  lower 
particle  velocities  than  those  assumed. 

As  in  the  previous  case,  the  current  standards  can  be  most 
effectively  compared  to  the  results  of  this  investigation  by  considering 
a  specific  example.  The  dimensions  and  rock  properties  used  in  the 
previous  comparison  will  be  assumed.  To  be  consistent  with  Equation 
5.10,  a  particle  velocity  of  9  in./sec  was  assumed  in  Equation  5.9.  The 
standoff  distance  required  according  to  Equation  5.10c  is  shown  in 
Figure  5.3  for  loading  densities  between  0.1  and  25  lb/ft3,  The  chamber 
separation  required  by  the  AESS  equation  Is  consistently  larger  than  chat 
required  by  Equation  5.9.  For  this  particular  example,  the  AESS 
standards  are  increasingly  conservative  with  loading  density.  The  AESS 
equation  appears  to  have  a  safety  factor  of  at  least  1.29  at  a  loading 
density  of  0.2  lb/ft3,  and  a  safety  factor  of  4.11  at  25  lb/ft3.  For 
loading  densities  of  2  to  3  lb/ft3,  the  AESS  equation  requires  standoff 
distances  approximately  twice  those  of  Equation  5.9. 

5.5  AIRBLAST  STANDOFF  DISTANCE 

To  determine  the  airblast  standoff  distance,  it  is  nccossary  to 
predict  accurately  the  pressure  at  the  tunnel  portal,  and  the  rate  of 
pressure  decay  with  distance  in  various  directions  beyond  the  portal. 
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Predictive  techniques  have  been  developed  to  determine  the  pressure 
generated  at  the  portal  for  elementary  geometries  if  the  loading 
densities  are  between  1  and  25  lb/ft3.  Typical  equations  were  given  in 
Chapter  1.  An  empirical  dimensionless  equation  which  delineates  the 
decay  of  the  peak  free- field  pressure  as  a  function  of  the  portal  exit 
pressure,  the  effective  diameter  of  the  tunnel,  and  distance  was 
developed  in  Chapter  4: 

P/PE  -  C^r/D.)-1*5  (5,11) 

where  P  is  the  peak  free-field  overpressure  (psi) 

PE  is  the  exit  pressure  at  the  tunnel  portal  (psi) 

D  is  the  effective  diameter  of  the  tunnel  (ft) 
and  Cx  is  a  constant  functionally  dependent  upon  direction. 

Based  upon  the  results  obtained  by  Gurke  and  Scheklinski 
(Reference  15),  Jenssen  (Reference  35),  and  this  investigation. 

Equation  5.11  will  predict  accurate  standoff  distances  near  and  along  the 
extonded  centerline  for  loading  densities  between  1  and  2.5  lb/ft3,  for 
exit  slopes  between  0  and  30  degrees.  At  lower  loading  densities,  the 
predicted  standoff  distances  will  be  increasingly  conservative.  Prudent 
use  of  Equation  5.11  is  required. 

Equation  5.11  is  based  upon  an  access  tunnel  which  exits  parallel  and 
tangent  to  smooth  terrain.  The  local  topography  can  significantly  alter 
the  standoff  distance  required.  In  general,  Use  attenuating  effects  of 
vegetation  and  the  natural  surface  roughness  tend  to  make  the  estimates 
conservative.  Palling  terrain  slopes  (with  respect  to  the  extended 
centerline  of  the  tunnel)  tend  to  decrease  the  pressure  levels. 
Conversely,  rising  slopes  may  tend  to  Increase  the  required  standoff 
distance.  Careful  consideration  should  be  given  to  any  surface 
characteristics  (hillsides,  valleys,  buildings,  etc.)  which  could  focus 
the  shock  wave. 

Equation  S. 11  may  not  be  valid  for  large  loading  densities.  As 
indicated  in  Figure  2.1,  at  TNT  loading  densities  greater  than 
0.07  lb/ ft3,  oxygen-deficient  combustion  products  result.  Secondary 
chemical  reactions  may  occur  in  the  access  tunnel  and/or  at  the  portal. 
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A  secondary  explosion  at  the  portal  would  increase  the  free* field  blast 
pressure  and  could  alter  the  directional  characteristics  of  the  peak 
overpressure.  The  magnitude  and  location  of  a  secondary  explosion  would 
depend  upon  the  chemical  composition  and  the  total  mass  of  explosive. 

Large  loading  densities  may  not  necessarily  result  in  large  blast 
pressures  at  the  portal.  Partial  or  complete  ground  shock  closure  of  the 
tunnel  exiting  the  storage  chamber  may  occur.  In  this  investigation,  a 
loading  density  of  25.3  lb/ft3  produced  partial  closure  and  reduced  the 
peak  free-field  blast  by  approximately  70  percent.  Tunnel  closure  has 
occurred  on  several  full*scale  tests,  and  is  sometimes  deliberately 
planned  to  prevent  venting.  As  a  result,  for  loading  densities  above  2.5 
lb/ft3,  Equations  2.1  and  5.11  should  be  used  with  caution  and  should  not 
be  employed  for  large  loading  densities. 

Existing  equations  cannot  be  used  to  determine  the  pressure  generated 
at  the  portal  of  complex  geometries,  because  these  equations  are 
empirical  relationships  which  define  the  peak  pressure  based  upon  total 
volume  of  the  complex  and  the  total  charge  weight.  As  a  result,  they  do 
not  explicitly  address  the  problem  of  shock  propagation  through  complex 
geometries.  If  the  blast  pressure  generated  at  the  portal  is  detorminod 
experimentally  for  a  complex  underground  geometry,  Equation  5.11  may  be 
used  to  predict  the  required  standoff  distance  if  the  loading  density  is 
not  sufficiently  large  to  produce  secondary  explosions  outside  the 
portal. 

The  airblast  standoff  distance  (D)  is  determined  in  the  AESS  by 
equations  of  the  following  form; 

D#  ~  C*  W*'3  <5.12) 

where  is  the  required  standoff  distance  along  a  line  which  is  at  a 
horizontal  angle  of  Q  degrees  from  the  extended  centerline 
and  W  is  the  not  explosive  weight  (lbs) 

is  a  constant  functionally  dependent  upon  the  direction  from  the 
portal. 
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This  equation  may  not  provide  reliable  estimates  of  the  required 
standoff  distance  in  all  cases,  because  it  applies  only  to  short, 
straight  tunnels  leading  directly  from  a  chamber,  and  does  not  assess  the 
geometrical  effects  on  shock  propagation  through  an  underground  complex. 

In  summary,  accurate  airblast  standoff  distance  can  be  determined  for 
elementary  tunnel/chamber  geometries  at  loa-'  .ug  densities  between  1  and 

2.5  lb/ft3.  At  high  loading  densities,  thr  current  AESS  equations 
predict  substantially  larger  standoff  distances  than  required.  For 
complex  geometries,  the  airblast  standoff  distance  can  be  determined  only 
if  experimental  data  for  that  geometry  is  available  to  predict  a  portal 
exit  pressure.  For  loading  densities  large  enough  to  produce  secondary 
explosions  inside  the  tunnel  complex  or  at  the  portal,  no  predictive 
techniques  are  currently  available. 

5 . 6  SUMMARY 

Equations  have  been  developed  in  this  study  to  predict  various  safety 
standoff  distances  required  as  a  result  of  explosions  in  underground 
magazines.  The  equations  for  ground  si  ck  address  magazine  geometrical 
effects,  explosion  decoupling,  and  the  properties  of  the  rock  medium. 
Additionally,  the  chamber  separation  distance  required  to  prevent 
symathetic  detonations  in  adjacent  chambers  by  rock  spalling  is 
functionally  dependent  on  the  impact  sensitivity  of  the  contents. 

The  standoff  distance  required  to  prevent  ground  shock  damage  to 
inhabited  buildings  is  related  to  the  peak  ground  shock  threshold 
velocity  of  Incipient  damage  to  the  structure.  The  equations  developed 
in  this  study  should  yield  conservative  standoff  distances  for  two 
reasons.  First,  they  are  based  on  upper  bound  strain  measurements.  In 
almost  all  instances,  the  actual  strain  field  will  be  less  than  the 
assumed  upper  bound  Second,  the  strain  data  was  obtained  around  a 
chamber  in  which  the  contents  were  detonated  essentially  instantaneously. 
In  full-scale  magazines,  sequential  detonation  of  the  contents  over  some 
finite  time  period  is  much  more  probable. 

The  current  AESS  equations  do  not  address  chamber  geometry  or 
decoupling  effects,  nor  do  they  have  a  capability  to  address  variations 
in  media  properties.  Essentially,  they  are  lumped  equations  which,  in 
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general,  will  predict  very  conservative  standoff  distances,  which  can 
result  in  increased  construction  costs  or  real  estate  requirements  (for 
buffer  zones) . 

The  data  obtained  in  this  study  corroborate  equations  proposed  by 
Skjeltorp  (Reference  4)  to  determine  the  standoff  distance  required  to 
prevent  damage  to  inhabited  buildings  by  airblast.  However,  these 
equations  are  valid  only  if  the  loading  density  is  sufficiently  small 
that  secondary  explosions  do  not  occur  near  or  outside  the  portal. 


Loading  Density  lb/ft 


Chamber  Separation,  ft 


Figure  5.1.  Calculated  storage  chamber  separation  required  to 
prevent  tunnel  closure  and  explosive  communication 
(by  spall),  compared  to  the  current  DOD  Safety 
Standard  (for  a  typical  design  in  granite). 
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Loading  Density  lb/ft 


Chamber  Separation,  ft 


Figure  5.2.  Calculated  storage  chamber  separation  required  to  prevent 
damage  to  contents,  compared  to  current  DOD  Safety 
Standards  (for  a  typical  design  in  granite). 
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Loading  Density  lb/ft 


Standoff  Distance,  ft 


Figure  5,3.  Calculated  standoff  distance  required  to  prevent 
damage  to  inhabited  buildings,  compared  to  current 
DOD  Safety  Standards  (for  a  typical  design  in  granite). 
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